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ABSTRACT 

In the framework of the recently proposed gluino-axion model, using the 
effective potential method and taking into account the top-stop as well as 
the bottom-sbottom effects, we discuss the CP-properties of the lightest 
Higgs boson, in particular its CP-odd composition, which can offer new 
opportunities at collider searches. It is found that although the CP-odd 
composition of the lightest Higgs increases slightly with the inclusion of the 
sbottom effects, it never exceeds %0.17 for all values of the renormalization 
scale Q ranging from top mass to TeV scale 







1 Introduction 



The radiative corrections to the masses of Higgs bosons in the MSSM have received much 
attention from the beginning, and an important step in the understanding of the MSSM Higgs 
sector was the significant modification of the tree level bound by the radiative corrections, 
dominated by the top quark and top squark loops 0, |||, ||]. The radiative corrections have 
been computed by using different approximations such as the effective potential [||, |||, ||], and 
diagrammatic |l], ||] methods. More complete treatment of these results include the complete one- 
loop on-shell renormalization ||, the renormalization group (RG) improvement for resumming 
the leading logarithms |@, ||, ||] , the iteration of the RG equations to two-loops with the use of the 
effective potential techniques fliol, [lj] , [l2|, |l^] and two loop on-shell renormalization [14, 15, |T(j |. 

In the recent literature, the studies on the radiatively induced CP violation effects and the 
theoretical predictions for the Higgs boson masses and couplings as functions of relevant minimal 
SUSY model(MSSM) parameters has been carried out in several directions. For example, in 
[17, 18], the implications of the presence of CP phases in the soft SUSY breaking sector allowing 
to the mixing of CP even and CP odd states were discussed. More recently the mass matrix of 
the neutral Higgs bosons of the MSSM has been calculated using the effective potential method, 
in the case of the small splittings between squark mass eigenstates in [jl9l , taking into account 
only the dominant top-stop contributions in [poj] and including the sbottom contributions in 



[21]. Additional contributions from the chargino, W and the charged Higgs exchange loops 
were computed in p2| . In |23fl , one-loop corrections to the mass matrix of the neutral Higgs 
bosons in the MSSM were calculated by using the effective potential method for an arbitrary 
splitting between squark masses, including the electroweak and gauge couplings and the leading 
two loop corrections. More complete treatment of the effective Higgs potential in the MSSM 
including the two-loop leading logarithms induced by top-bottom Yukawa couplings as well as 



those associated with QCD corrections by means of RG methods were performed in [24|. 

It is a well-known fact that in the supersymmetric (SUSY) extensions of S.M, apart from 
the physical phases 8c km and 9qcd, also existing in the SM, there appear novel sources of 
CP violation via the phases of the soft supersymmetry breaking mass terms J25|. Besides their 
contribution to the known CP-violating observables such as electric dipole moments of the 
particles ^(| |?]] , these phases also induce CP violation in the Higgs sector 0, |l8| [H^, |2^, |||, |t| . 
Moreover, the SUSY theories which are designed to solve the hierarchy problems, possess two 
hierarchy problems: One concerning the strong CP problem, also existing in the S.M, whose 
source is the neutron EDM exceeding the present bounds by nine orders of magnitude |3(|] and 
the other is the \x puzzle, concerning the Higgsino Dirac mass parameter /j, which follows from 
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the superpotential of the model. A simultaneous solution to these two hierarchy problems can 



be achieved by the gluino-axion model |3l|, 32, 33 1 with a new kind of axion [34, 35, 36, |37|, which 



couples to the gluino rather than to quarks. Besides, the low energy theory is identical to the 
MSSM with all sources of the soft SUSY phases. Due to all these abilities of the model, in the 
analysis below we shall adopt its parameter space. 

In the recent literature, the radiative corrections to the Higgs masses and mixings, dominated 



by the top-stop contributions have been studied in gluino-axion model [39|. In Ref. |39|| , which 



was based on Ref. [20], the mass squared matrix of Higgs scalars involves the one loop function 
g(x, y), as well as the other scale-independent terms. The function g(x, y) is related to f(x, y) in 
such a way that g(x, y) = f(x, y) — log^, where f(x, y) is a scale-dependent loop function whose 
expression is given by f(x,y) = —2 + log^ + ~^log&. However, since the explicit dependence 
of Q in the one- loop function f(x,y) is actually cancelled by the explicit Q dependence of the 
function g(x,y), then, g(x,y) = —2 + ^jz^log^ does not have an explicit dependence on the 
renormalization scale Q, unlike f(x,y). Therefore, in both of the works of [20, |3^] the elements 
of the mass squared matrix of Higgs scalars, which depend only on the scale-independent function 
g(x, y) as well as the other scale-independent terms, do not have an explicit dependence on Q. 



Actually, in [20, |39| the one- loop bottom-sbottom contributions to AV are not included, and the 
terms proportional to fi, At are obtained by neglecting the D-terms in the stop masses, to gain 
independence of the renormalization scale Q, since the D-term contributions to the squark masses 



are quite small. On the other hand, in Ref. [21] the one-loop bottom sbottom contributions as 
well as the top-stop contributions to radiative corrections are taken into consideration. Thus, 
the elements of the mass squared mass matrix of the Higgs scalars not only include the scale- 
independent function g(x, y) with various parameters, but have an explicit dependence on the 
renormalization scale through the scale dependent function f(x,y), and the scale dependent 
logarithmic terms, stemming from the additional contributions to the radiative corrections. 
The main purpose of this work is to analyze the CP violation effects on the lightest Higgs 



boson in the framework of the gluino-axion model using the recent experimental data [38], by 
the inclusion of the bottom-sbottom effects, as well as the top- stop contributions. As the (In- 
dependence is taken into consideration, we particularly adress the issues, whether or not the 
various renormalization scales, all being around the weak scale, would lead us to a large amount 
of CP violation opportunities, and, investigate whether one can find an appropriate limit of 



reasonable agreement with the scale-independent results [39|. We will base our calculations to 



those of Ref. [|2lf| , However, we differ from [21] in the sense that in our analysis, all the chosen 
parameters are specific to the gluino-axion model, namely all the soft mass parameters in this 
theory are fixed in terms of the [i parameter. 
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The organization of this work is as follows: In Sec. 2, starting from the Higgs sector structure 
of the gluino-axion model, we compute the (3 x 3) dimensional mass matrix of the Higgs scalars 
whose elements are expressed in terms of the parameters of the model under concern. In Sec. 3, 
we make the numerical analysis to study the CP violation effects on the lightest Higgs boson, 
and analyze the influence of Q on its CP-odd composion, as well as on fi. We conclude in Sec. 4. 

2 The model 

In the gluino-axion model, the invariance of the supersymmetric Lagrangian and the invariance 
of all supersymmetry breaking terms under U(1)r are guaranteed by promoting the ordinary 
ji operator to a composite operator containing the singlet composite superfield S with unit R 



charge [31]. When the scalar component of the singlet develops a vacuum expectation value 
around the Peccei- Quinn scale v s ~ 10 11 GeV, an effective ^-parameter ~ a TeV is induced 
such that 

H = v 2 s /Mpi x e - id Q CD/3 ~ a TeV x e - ie ^ CD/3 (1) 

where Oqcd is the effective QCD vacuum angle. Therefore, the vacuum expectation value of the 
singlet serves for two important purposes for the model under concern: Its magnitude determines 
the scale of supersymmetry breaking and its phase solves the strong CP-problem. 



The effective Lagrangian at low-energy is given by [31]: 

£mssm = & m qQ + u- ] M 2 uC u- + d^M 2 J- + Vm 2 L + e^M 2 c e- 

+ [A U Q ■ H u u c + A d Q ■ H d d c + A e L ■ H d e c ] + h.c.} 

+ M 2 H JH U \ 2 + M 2 Hd \H d \ 2 + fa BH U ■ H d + h.c.) 

+ {M 3 A^A^ + M 2 A l 2 A* 2 + M 1 AiAi + /i.c.}, (2) 

The soft terms of the low energy Lagrangian in the gluino-axion model are identical to those 
in the general MSSM except for the fact that the soft masses are all expressed in terms of the 
H parameter through appropriate flavour matrices except for the fact that the soft masses are 
all expressed in terms of the n parameter through appropriate flavour matrices. The flavour 
matrices form the sources of CP violation and intergenerational mixings in the squark sector. 
The phases of the trilinear couplings (At,d,e)> the gaugino masses (Ms ; 2,i), and the effective 
/^-parameter defined in (||) are the only phases which can generate CP violation observables. 

It is known that the dominant contributions to the one-loop radiative corrections come 
from the the top quark and top squark loops as long as tan/3 < 50, as in the CP-conserving 
case JI], ||, ||, H]. However, for a more sensitive calculation, we will take into account of the 
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contributions from the bottom- sbottom, as well as the top-stop quark loops. For convenience, 
we set the soft SUSY breaking scalar-quark masses as Mq = = Mj, and the squark trilinear 
couplings as A t = Ay,. Then, the explicit expressions for the mass parameters in @) defined as 
follows: The top (bottom) squark soft masses are given by: 

M| = k 2 Q M 2 , (3) 

where kQ is a real parameter. The top(bottom) squark trilinear couplings read as: 

A t = fxkl (4) 

where kt is a complex parameter. The tree level Higgs soft masses are defined by: 

M 2 Hu = y u |/i 2 | , M 2 Hd = y d \ix 2 \ , fi B = H 2 (^ + , (5) 

D Q 



where m 2 ~ v 2 is a natural choice as discussed in [31|. Here y u and are real parameters, and 
k^ is a complex parameter determining the phase of the B parameter, which can be identified 
with the relative phase of the Higgs doublets p0|. 



















f Hi 









After electroweak breaking the Higgs doublets in (^) can be expanded as 
H d = 

Hu = [ h} ) = 7! [ , u + 2 u +^ 2 ) ■ (6) 

where tan/3 = v u /v d as usual, and the angle parameter 6 is the misalignment between the two 
Higgs doublets. 

We follow the effective potential method for computing the one-loop corrected Higgs masses 
and mixings. As usual, the entries of the Higgs masses and their mixings can be calculated up to 
one loop accuracy by the second derivatives of the effective potential with respect to the Higgs 
fields. 

( d 2 V \ 

M 2 = , where Xi S 5 = {0i, (/) 2 ,^i, 9? 2 } , (7) 

\dXidXjJ Q 



where V = Vq + Vi-i oop is the radiatively corrected Higgs potential [21], and as we mentioned 
before we take into account the top-stop and bottom-sbottom loop corrections. The stop and 
sbottom mass-squared eigenvalues are given by: 

m L = - A M 2 z C2 P + m 2 + k 2 Q \^\ 2 ^^l , (8) 
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m? i 2 = ~M 2 z C2i3 + m 2 b + k 2 Q \ f i\ 2 T^l , (9) 

with, 

Af = J (bfy - ^M 2 fc% + m 2 M 2 (\h\ 2 + tf - 2\k t \t^ c m ) , (10) 



where 0/3=003/3, c Vkt =cosipkt, sp=sw./3, s,p kt =smipk t , tg=tan/3, t^ 1 =cot/3. The stop and sbot- 
tom mass splittings depend explicitely on the total CP violation angle ipkt such that 

ip kt = Arg[M f *] = Arg[A: t ] , (12) 

where kt has been defined in (Q). 

The (3x3) dimensional Higgs mass-squared matrix can be expressed as: 

Mfc 2 + M 2 s 2 + A u -(M 2 +M]) S/3C/3 + A 12 A 13 
-(Mf + Mi) spcp + Ai 2 Mfs 2 + M 2 c 2 + A 22 A 23 I (13) 



Ai 3 A 23 Mjf + A 33 

in the basis B = {0i, 2 , sin^i + cos/3^} using (||). 

The elements of the radiatively corrected mass-squared matrix read as: 
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where = 8n $ v 2 an d Q is the renormalization scale in the MS scheme. In the above expressions, 



ftn = p (\hK kt +t f3 1 ) n t2 = /i(|fct| + t j 9 1 <v w ) 

= /i(|fet|<V fct H- *j3 ) =/i(|fct| + */9<V w ) 
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and, 



X t = Z t (M w , M z )c 2f 3 , X b = -Z b (M w , M z )c 2p . (27) 

where 

Z t (M w ,M z ) = (Im^-^mI) 2 , 

Z b (M w ,M z ) = (^M^-^M 2 ) 2 , (28) 

and the functions f(m 2 .- , ,m 2 .- .) and q(m 2 - ,,m 2 ,- ,) are given by: 
J y tl(6i)' t 2 (b 2 )> yv tl(6i)' t 2 (b 2 )' fa J 

m~ - m~ - 

f(m 2 ~ m 2 - ) - -2 + log l( l} 2( 2) + 

m 2 , 7 , + m 2 . m 2 , r , 

-MM ^lliog^M (29 ) 

m- ,- . — m~ ,r , m- ,- , 
*2(&2) ti(!n) t 2 (b 2 ) 

TTl 2 ~ ~\~ TTl 2 - 7Yt 2 ~ 

I 2 2 \ o i *a( 6 a) *i(*>i) , fo(&a) / on \ 

g{m~ r m ; ( - t = -2 H o-^^ log 9 (30) 

ti(6i)' <2(fe) / 2 2 & 2 

ta(6a) *i(6i) ti(bi) 
We diagonalize the Higgs mass-squared matrix (12) by the similarity transformation 

UM 2 1Z T = diag(m^ , m\ 2 , m 2 h:j ) , (31) 

where VSRF = 1, and we define /13 to be the lightest of all three. 

One of the most important quantities in our analyses is the percentage CP composition of 
a given mass-eigenstate Higgs boson. The percentage CP compositions of the Higgs bosons in 
terms of the basis elements are defined by 



Pi 



100 x \K U \ 2 ; i = 1,2,3. (32) 



where pi, pi and p^ correspond respectively the <f>i, 02, sin/3</?i + cos 0cf2 components of the 
Higgs boson under concern. In the following, one of the main concerns will be the CP-odd 
composition p$, of the lightest Higgs boson as it can offer new opportunities at colliders for 



observing the Higgs boson [29, 40, 41, 42 1. We will discuss the dependence of p% on different 



renormalization scales, as well as the interdependence of Q-p, using the the previous theoretical 



, as well as the experimental results [38]. 
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3 Numerical Analysis 



In this section, we mainly adress the phenomenological consequences of explicit radiative CP 
violation effects in the lightest Higgs boson, whose CP-odd composition as well as its mass are 
of prime importance in direct Higgs boson searches at high-energy colliders. As a result of the 
standart model Higgs boson searches at LEP, the lower bound on the lightest Higgs mass is 
115 GeV (and correspondingly tan/3 > 3.5) [f&jjl . On the other hand, theoretically the lightest 



Higgs boson mass can not exceed 130 GeV for large tan/3 [16]. Therefore, from the searches 
at LEP2, the lower limit on mass of the SM Higgs boson excludes the substantial part of the 
MSSM parameter space particularly (for m t = 175 GeV), at small tan/3 (tan/3 < 3.5) [p8| . 

Being a reflecting property of the underlying model, all the soft masses are expressed in 
terms of the p parameter, and since the p parameter is already stabilized to the weak scale as a 
consequence of the naturalness, all dimensionless quantities are expected to be 0(1). Therefore, 
we take all the dimensionless quantities in the 0(1). In our numerical analysis, we study two 
specific values tan/3, namely tan/3 = 4 and tan/3 = 30, to analyze the CP violation effects on 
the lightest Higgs boson in the the low and high tan/3 regimes. 



Q_ < " 1 




Q.P> 




Figure 1: The CP-odd composition (p^) of the lightest Higgs boson as a function of (fkt, for 
tan/3 = 4 (left panel), and tan/3 = 30 (right panel), where Q varies from 175 GeV to 1200 GeV. 



Shown in Fig. 1, is the ifkt dependence of ps, designating the percentage CP-odd composi- 
tions of the lightest Higgs boson (/13), for tan/3 = 4 (left panel), and tan/3 = 30 (right panel). 
In both panels Q changes from 175 GeV to 1200 GeV, and \p\ from 100 GeV to 1000 GeV in 
the full (fkt range. Therefore, the variation of p^ is represented by many points in the parameter 
space which correspond to different values of Q, and different values of p. As is seen from both 
panels of Fig. 1, the maximal value of p% occurs at %0.025 for tan/3 = 4, and at %0.17 for 
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tan (3 = 30. Moreover, a comparative look at both windows shows that, although P3 increases 
relatively with the increasing tan/3, it never exceeds %0.17 in the full ifkt range, for all values of 
Q changing from from 175 GeV to 1200 GeV. Compared to its CP-even compositions, which 
form the remaining percentage, this CP-odd component is extremely small to cause observable 
effects. It may, however, be still important when the radiative corrections to gauge and Higgs 
boson vertices are included 

The analysis of Fig. 1, gives a general idea of the variation of ^3 in the full (p^ range. 
However, since the parameter space is quite large, for a better understanding of the properties 
of the model under concern, we focus on the different portions of parameter space corresponding 
to the different values of Q in Fig. 2, Therefore, in the left panel of Fig. 2, we choose three 
specific values of Q corresponding to 325 GeV ("+"), 600 GeV ("x"), 1000 GeV ("."), for 
tan (3 = 4, and analyze the variation of p% in the full (p^t range, for all values of [i changing from 
100 GeV to 1000 GeV. On the other hand, we carry out the same analysis for tan/3 = 30 in the 
right panel of Fig. 2. However, for convenience, starting from a lower value of Q, we let 190 GeV 
("o"), 325 GeV ("+") , 1000 GeV (".") values of Q for this case. 
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Figure 2: The CP-odd composition (/%) of the lightest Higgs boson as a function of ifkt, 
for Q = 325 GeV ("+"), 600 GeV ("x"), 1000 GeV ("."), for tan/3 = 4 (left panel), and 
Q = 190 GeV (V), 325 GeV ("+"), 1000 GeV ("."), for tan/3 = 30 (right panel). 



In Fig. 2, those portions of the parameter space which belong to different values of Q are 
presented by curves ("+", "x", "o" ), and those curves correspond to different values of 
fi. A comparative look at Q = 325GeV ("+"); and Q = 600 GeV ("x") scales, in the low tan/3 
regime (left panel), shows that when Q = 325GeV ("+") and tan/3 = 4 (left panel), there are 
very few curves, since the lower bound of \x (corresponding to Q = 325GeV) starts from a rather 
higher value of n (/i = 800GeV), and, it is not possible to find any allowed region below this 
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value (p < 800GeV), since the parameter space is constrained by the existing LEP bound on 



the lightest Higgs boson mass [38]. Moreover, even at this value of p, one can not find solutions 
in the full (fkt range, as will be shown in Fig. 3. Therefore, the allowed range of the parameter 
space is quite restricted for this case. At Q = 600GeV and tan/? = 4 (left panel), p3 decreases 
relatively as compared to the former case. However, as the allowed range of p gradually gets 
widened, the number of curves increase (for instance, the lower bound of p is 750GeV for this 
scale). On the other hand, p^ is the quite small at Q = lOOOGeV, and tan/? = 4 (left panel), 
as compared to the other two scales (Q = 325GeV, Q = 600GeV). However, there is a gradual 
enlargement in the allowed range of p (The allowed range of p for all values of Q will be discussed 
in detail, in analyzing the Q — p interdependence). The variation of p% — (fkt is similar in the 
high tan/3 regime (right panel of Fig. 2). Considering the portion of the parameter space which 
corresponds to Q = 190GeV ("o"), one notes that the parameter space is quite widened, even 
as compared to a larger scale, Q = 325GeV ("+"), in the low tan/3 regime. On the other hand, 
a comparision between the left and right panels, at Q = 325GeV, also shows that, the number 
of curves, namely the allowed range of p increases gradually, as well as p$, in the high tan/3 
regime (right panel). 

From the analysis of Figs. 1, 2, we can conclude that the allowed range of p, as well as p% 
gradually increases with the increasing tan /?. However, as the p% — p dependence is taken into 
consideration, it is seen that some differences arise in connection with the different values of Q. 
Therefore, to understand the properties of different values of Q, and in particular their influences 
on p3, as well as on p, we carry out the same analysis from a slightly different perspective in 
Figs. 3, 4, 5 focusing on the behaviour of the small, moderate and large values of Q, under 
different values of p, in both regimes 

Shown in Fig. 3, is the variation of the CP-odd composition (^3) of the lightest Higgs 
boson as a function of cpkt, for tan/? = 4 (left panel), and tan/3 = 30 (right panel), when 
p = 800 GeV. The three different curves present three different values of Q corresponding to 
325 GeV (top curve) , 600 GeV (middle curve) , 1000 GeV (bottom curve), respectively. The same 
analysis has been carried out for p = 900 GeV in Fig. 4, and, for p = 1000 GeV in Fig. 5. As is 
seen from the left panel of Fig. 3, at Q = 325 GeV, and p = 800 GeV (the lower bound of p for 
this scale), there are solutions only in a small portion of parameter space in the full (fkt range, 
which corresponds to [4-7r/5, 6tt/5] interval, since one can not find any allowed region satisfying 
the experimental constraint |38| beyond this interval. At Q = 325 GeV, and p = 900 GeV (left 
panel of Fig. 4), the allowed range of ipkt slightly widens. One notes that, for Q = 325 GeV, the 
experimental constraint is satisfied in the full (fkt range, only at p = 1000 GeV (left panel of 
Fig. 5). On the other hand, as the renormalization scale increases, for instance, at Q = 600 GeV, 
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Figure 3: The CP-odd composition (/J3) of the lightest Higgs boson as a function of (fkt, for 
tan/3 = 4 (left panel), and tan/3 = 30 (right panel), and p = 800 GeV, where the top, middle and 
bottom curves correspond to Q = 325 GeV ("+"), 600 GeV ("x"), 1000 GeV ("*") respectively. 

the allowed range of ifkt gradually gets widened, and at Q = 1000 GeV, there are solutions in 
the full ifkt range. Similar observations can be made for the tan/3 = 30 case. However, as is 
seen from the right panels of Figs. 3, 4, 5 that the experimental constraint is satisfied for all 
values of (fkt in this regime. 
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Figure 4: The same as Fig. 3, but for p = 900 GeV 



Moreover, as is shown in the left panel of Fig. 4, for Q = 325 GeV, the maximal value of 
P3 occurs at ~ 0.004 for p = 800 GeV. Then, increasing slightly, and it reaches to ~ 0.016 at 
[i = 900 GeV (the left panel of Fig. 5). Finally, the maximum value of P3 occurs at ~ 0.022 at 
fi = 1000 GeV (left panel of Fig. 6 ). On the other hand, p% slightly decreases for Q = 600 GeV, 
as compared to the former case (Q = 325 GeV), and like the former case, its p^-p dependence is 
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Figure 5: The same as Fig. 3, but for p = 1000 GeV 



of the same kind (the larger p, the larger ps). In passing to Q = 1000 GeV, one notes from the 
left panel of Fig. 3 that, P3 decreases relatively as compared to Q = 325 GeV, and Q = 600 GeV 
cases, and it occurs maximally at ~ 0.0005 at p = 800 GeV. However, unlike the Q = 325 GeV 
and Q = 600 GeV cases, ps (corresponding to Q = 1000 GeV) is smaller at p = 900 GeV than 
that of p = 800 GeV, and in fact, it reaches to its minimal value at p = 1000 GeV (which will 
be seen in detail in analyzing the ps — p interdependence of Q = 1000 GeV). 

Therefore, the comparative analysis of Figs. 3, 4 and 5 shows that, there are differences in 
the P3~p dependence of Q in both low and high tan (3 regimes, in the sense that, although p% 
increases with increasing p for Q = 325 GeV, and, Q = 600 GeV, its dependence on p changes 
in passing to the higher renormalization scales (for instance at Q = 1000 GeV, p% decreases with 
increasing p). Therefore, one can deduce from the above analysis that, there should be critical 
values of Q in both low and high tan /3 regimes, beyond which the pz-p dependence changes. To 
determine these portions of the parameter space, it is necessary to discuss the variation of p^, 
with Q, which will be carried out in the following. 

Before analyzing the dependence of Q on p^ in the full ifkt range from a general point of 
view, we focus on a smaller portion of a parameter space, for clearance. Therefore, in Fig. 6, 
we choose a particular value of ipkt, namely, we set (fkt = 37r/2. Moreover, as Q changing from 
175 GeV to 1200 GeV, we let 975GeV ("*"), 875GeV ("x"), 775GeV ("+"), 675GeV ("o") 
values of p, for tan/3 = 4 (left panel), and tan/3 = 30 (right panel). One notes that, there is no 
allowed region in the parameter space corresponding to p = 675GeV for tan/? = 4 (left panel). 
In fact, it is not possible to find any region in the parameter space below p < 700 for tan/3 = 4 
(this bound is p < 500 for tan (3 = 30), due to the experimental constraint on the lower bound of 



the lightest Higgs mass [38]. Thus, in choosing the different values of p, we try to form the most 
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Figure 6: The CP-odd composition (p^) of the lightest Higgs boson as a function of the renor- 
malization scale Q, at ipk t = 37r/2, for tan/3 = 4 (left panel), and tan/3 = 30 (right panel), 
with different values of p corresponding to p = 975GeV ("*"), 875GeV ("x"), 775GeV ("+"), 
675GeV ("o") 



suitable combination. In Fig. 6, for the ease of following we cut the vertical axes p% = 0.02% for 
tan/3 = 4 (left panel), and p$ = 0.06% for tan/3 = 30 (right panel). However, as will be seen 
in Fig. 7, it actually extends to 0.025% for tan/3 = 4, and 0.18% for tan/3 = 30 in the full <pkt 
interval. Although the analysis has been carried out for a particular value of (fkt {(fikt = Stt/2), 
one notes that there are two portions in the parameter space. In the first portion, p% decreases 
with decreasing p, then after a certain value of Q, in the second portion, p% starts to increase 
slightly with decreasing p. For instance, as is seen from the left panel of Fig. 6, in the first 
portion of the parameter space (Q ^ 925GeV), the maximal value of p% occurs at p = 975GeV 
("*"), and it decreases slightly with the decrease in p. On the other hand, in the second portion 
of the parameter space (Q j£ 925GeV), the maximal value of p% is smaller than that of the first 
region. However, p% is maximal at p = 775GeV ("+"), and, unlike the first portion of the 
parameter space, p% lessens as p increases. For instance, p% is minimal at p = 975GeV ("*"), in 
the second portion. Similar observations can be made also for the high tan/3 regime. Clearly, 
with the analysis of Fig. 5, we focus on a particular value of (fkt, and try to understand the 
P3~p dependence of Q, considering different values of p. In the following, we will carry out the 
same analysis from a more general point of view, which will provide us to determine the critical 
values of Q more precisely, and to understand the properties of the model under concern in more 
detail. 

In Fig. 7, we show the variation of p% with Q for tan (3 = 4 (left panel) , and tan (3 = 30 
(right panel) in the full (p^t range. In each panel, Q ranges from 175 GeV to 1200 GeV, p 
from 100 GeV to 1000 GeV. Here, the vertical lines represent different values of ip^t, as well 
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Figure 7: The CP-odd composition (p^) of the lightest Higgs boson as a function of the renor- 
malization scale Q, for tan/3 = 4 (left panel), and tan/3 = 30 (right panel). 

as different values of p, as was indicated in Fig. 6. As is seen from the left panel of Fig. 7, 
the maximum value of p% (~ 0.025%) occurs at Q = 275 GeV. For larger values of Q, 
decreases gradually until Q = 925GeV, as the super symmetric spectrum decouples. One notes 
that, beyond this point p% starts to increase slightly. For smaller values of Q, however, the 



parameter space is constrained by the existing LEP bound on the lightest Higgs mass [38]. That 
is, p3 gets smaller until Q = 175 GeV, and it is not possible to find any region in the parameter 
space below this value (Q ^ 175 GeV), for tan/3 = 4, since this region is completely disallowed 
by the experimental bound |3Sj| . On the other hand, as is shown in the right panel of Fig. 7, 
the maximum value of p% (~ 0.17%) occurs at Q = 175 GeV for tan/? = 30. Similar to the low 
tan (3 regime, p% decreases with increasing Q until Q = 775GeV, and beyond this point there is 
a small increase in p% in the admitted range of Q. 

Following the analysis of Figs. 3-7, we would like to emphasize that in the first portion of the 
parameter space (175 ^ Q ^ 925 for tan/3 = 4, and 175 < Q ^ 775 for tan/3 = 30), as p increases 
with increasing p% (Figs. 3-6), p3 decreases with increasing Q (Figs. 6-7). On the other hand, in 
the second portion of the parameter space (925 ^ Q ^ 1200 for tan (3 = 4, and 775 ^ Q ^ 1200 
for tan/3 = 30), p3 gradually increases with increasing Q (Figs. 6-7), and since this increase 
in /?3 is compensated by the decrease in p, unlike the first region, p decreases with increasing 
P3 (see Fig. 3-6). Therefore, the interdependence of p — p% changes in this region (see Fig. 6). 
However, the maximal values of p% never exceeds ~ 0.005% for tan/3 = 4, and ~ 0.02% for 
tan P = 30, in the second portion of the parameter space. Obviously if the lightest Higgs boson 
were carrying large enough CP-odd composition, this would bring about new opportunities for 
observing the lightest Higgs in the near future. However, one notes that as the CP violation 
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effects are induced only radiatively in the MSSM and in the model at hand, it is clear that any 
would-be observation of a large CP violating composition of the lightest Higgs boson at LHC, 
or NLC, or TESLA excludes the MSSM in general, and the model under concern in particular. 
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Figure 8: The dependence of the renormalization scale Q on |/i| for tan/3 = 4 (left panel), and 
tan /3 = 30 (right panel) . 



In Fig. 8, we show the interdependence of Q on \fj>\ for tan f3 = 4 (left panel), and tan/3 = 30 
(right panel). In each panel Q varies from 175 GeV to 1200 GeV, and /i from 100 GeV to 
1000 GeV. As is seen from the left panel, in the first portion of the parameter space, starting 
from Q = 175GeV, the range of Q gets widened up to Q = 925GeV, whereas the lower allowed 
bound of /j, changes from 850 GeV to 700 GeV in this interval. For instance, at Q = 325 GeV 
and tan/3 = 4, the lower bound of /i starts from 800 GeV (see Fig. 3), and, at 600 GeV, it 
decreases to 750 GeV in the same regime, as will be shown in Fig. 9. On the other hand, in 
the second portion of the parameter space (925 < Q < 1200), the allowed range of /i remains 
constant. That is, all values of /j, changing from 700 GeV to 1000 GeV are allowed for Q > 925 
(for instance, the lower bound of /i, corresponding to Q = lOOOGeV starts from 700 GeV, as will 
be indicated in Fig. 10). Q-/U interdependence is similar at large tan/3 (right panel). However, 
in this regime, as the range of Q gradually widens in the 175 < Q < 775 GeV interval, the lower 
bound on /j, changes from 600 GeV to 500 GeV. Above this interval {Q > 775), all values of /u, 
from 500 GeV to 1000 GeV are allowed. For instance, the lower bound of (i, corresponding to 
600 GeV starts from /j, = 550 GeV at tan /3 = 30, whereas it is n = 500 GeV, for Q = lOOOGeV 
in the same regime, beyond which all values of fi are allowed (see Figs. 9-10). 

We would like to remind that in the present analysis, the allowed range of the parameter 
space is obtained by imposing the recent LEP constraint on the lightest Higgs mass (mh 3 ) which 



requires m^ 3 > 115GeV [38|. As is seen from both panels of Fig. 8, the lower bound on /i 
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starts from 700GeV for tan/3 = 4 (left panel), whereas this bound decreases to 500GeV for 
tan/3 = 30 (right panel), and it is not possible to find any region in the parameter space below 
these values (/i < 700GeV for tan/3 = 4, and p < 500GeV for tan/3 = 30), since these regions 
of the parameter space are completely disallowed by the existing LEP bound on the lightest 
Higgs mass j38|. Moreover, it is worthwhile of mentioning that the present analysis includes the 
dominant top-stop contributions, as well as the bottom-sbottom effects, whereas the two-loop 
corrections have not been taken into account. In this case, the approximation used here for the 
calculation still contains a theoretical error of several GeV |43[[ . Therefore, the allowed range 
might change mildly, in case of a higher precisional calculation. 

With the former analysis, we have studied the properties of the model under concern, as 
the Q-dependence is taken into the consideration. At this point, we come back to the question 
mentioned in the purpose of this work (whether one can find an appropriate limit of reasonable 
agreement with the scale independent results |^]]). Therefore, in Fig. 9, we focus on two 
particular values of Q, namely Q = 600 GeV, and Q = 1000 GeV, and investigate their p% — p 
dependence in both high and low tan/3 regimes. 
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Figure 9: The CP-odd composition (p^) of the lightest Higgs boson as a function of \p\ for 
Q = 600 GeV, when tan/3 = 4 (left panel), and tan/3 = 30 (right panel) 



Depicted in Fig. 9 is the \p\ dependence of p3 for tan/3 = 4 (left panel), and tan/3 = 30 (right 
panel) at Q = 600 GeV. In Fig. 10, the same dependence is shown for a larger scale, namely for 
Q = 1000 GeV. In both panels, the vertical lines correspond to different phases in the full ipkt 
range (see Figs. 3, 4, 5). A comparative glance at both panels shows that, in both high and low 
tan/3 regimes, the P3 — p dependence at Q = 600 GeV scale differs from that of Q = 1000 GeV. 
In fact, the ps — ifkt dependence corresponding to these scales has been analyzed in Figs. 3, 4, 5, 
and the critical values of Q (beyond which p% — p behaviour changes) has been determined in 
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Figure 10: The same as Fig. 9 , but for Q = 1000 GeV. 



Figs. 6, 7. Therefore, such kind of behaviour is expected, remembering that in the first portion 
of the parameter space (175 < Q ^ 925 for tan/? = 4, and 175 < Q ^ 775 for tan/3 = 30), 
p increases with increasing p% (Fig. 7), whereas it decreases with increasing p% in the second 
portion of the total Q range (925 <Q< 1200 for tan (3 = 4, and 775 < Q < 1200 for tan f3 = 30). 
As is seen from the left panel of Fig. 9, at Q = 600 GeV, and tan (3 = 4, the maximum value of 
P3 (~ 0.005%) occurs at fj, = 1000 GeV, and for smaller values of p, p^ decreases gradually. In 
passing to tan/3 = 30 case (right panel), one observes that as p3 increases slightly, the allowed 
range of p widens (the lower bound on p decreases to p = 550GeV). On the other hand, as is 
seen from the left panel of Fig. 10, at Q = 1000 GeV, and tan/3 = 4, the maximal value of p% 
(~ 0.001%) occurs at p = 700 GeV, and unlike the Q = 600 GeV case, for larger values of p, p^ 
decreases gradually. The variation of p% with p, for tan (3 = 30 (right panel of Fig. 10) is similar 
to what we have observed in the tan/3 = 4 regime (left panel of Fig. 10). However, there is a 
gradual increase in p%, as well as the allowed range of p, for tan/3 = 30. 

We would like to note that the variation of p3 with p (Fig. 10), which corresponds to 
Q = 1000 GeV, is very similar to what we have found in j39|. This particular result shows that 
the scale dependence of the radiative corrections is sufficiently suppressed for the renormalization 
scale Q = 1000 GeV. On the other hand, consideration of various renormalization scales, all 
being around the weak scale, lead us to a wealth of CP violation opportunities. This is especially 
suggested by Fig. 9, where the percentage CP-odd composition of the lightest Higgs gradually 
increases with p. 

As mentioned above, among the values of Q changing from 175 GeV to 1200 GeV, of partic- 
ular interest is the Q = lOOOGeV where the comparision between the results of this work and 
that of [E^] is possible, and a reasonable agreement with (3J3] can be found. Shown in Fig. 11, 
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Figure 11: The dependence of the CP-odd composition (^3) of the lightest Higgs boson on tp^, 
corresponding to Q = lOOOGeV, for tan/3 = 4 (left panel), and tan/3 = 30 (right panel) 



is the tfkt dependence of the percentage CP-odd compositions of the lightest Higgs boson (^3), 
corresponding to Q = 1000 GeV, for tan/3 = 4 (left panel), and tan/? = 30 (right panel). As 
both panels of the figure suggests the maximal value of p% occurs at ~ 0.001% for tan (3 = 4 (left 
panel), and increasing slightly it reaches to ~ 0.009% for tan/? = 30 (right panel), in the full 
(fkt range. As has been shown in Figs. 3,4,5 each curve in the figure represents different values 
of (a. For instance, in the tan (3 = 4 regime (left panel) the curve on the top corresponds to the 
lower bound of p (p = 700GeV), for which p% reaches its maximal (see Fig. 10), and similarly 
from the top curve to the bottom with the increase in p,, p% decreases. Similar to observations 
made for the left panel, one can discuss the high tan /? regime (right panel) , where the top curve 
corresponds to the lower bound of p (p = 500GeV), for tan/? = 30 (see Fig. 10). One notes 
that, the p% reaches maximally to ~ 0.009% in the full ifkt range, for all values of tan (3 changing 
from 4 to 30. 

As is explained in the Introduction, considering only the dominant top-stop quark loops, 
and being the results of Q-independent, the variation of P3 with ip^t is analyzed in p9| , and 
it has been found that the maximal value of ps never exceeds 0.0013% for all values of tan/? 
changing 4 to 30 in the full <pkt range (the maximal values of P3 ~ 0.0003% for tan (3 = 4:, 
and ps ~ 0.0013% for tan/3 = 30 in |39[]) The dependence of ps on ipkt (see Fig. 10), which 
corresponds to Q = 1000 GeV, is very similar to that of [39]. However, one notes that maximal 
value of P3 slightly increases in this work, as compared to the former analysis. 

Among the particles contributing to the one-loop radiative corrections, the dominant ones 
come from the top quark and top squark loops, provided that tan (3 < 50, (for which case the 
bottom Yukawa coupling is quite small to give significant contributions). The Yukawa interac- 
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tions due to scalar bottom quarks can be significant only for very large tan (3 values. Therefore, 
by the inclusion of the bottom-sbottom loops p% increases slightly as expected, however, it never 
exceeds 0.009%, for all values of tan (3 changing from 4 to 30. 



f) m 




Figure 12: The dependence of the mass of the lightest Higgs boson (m^ 3 ) on (p^ t , corresponding 
to Q = lOOOGeV, for tan/3 = 4(left panel), and tan/3 = 30 (right panel) 

Finally, depicted in Fig. 12, is the variation of /?3 with m/ l3 , corresponding to Q = 1000 GeV, 
for tan /3 = 4 (left panel) , and for tan (3 = 30 (right panel) . The decreasing curves in each panel 
represent different values of ipkt, and p as well. The maximal value of pz occurs at ~ 0.001%, 
and correspondingly at p, = 700GeV (see Figs. 10, 11), Then, it decreases rapidly with the 
increasing mass, and decreasing p. For instance, it reaches to ps ~ 0.0005% at p = 800GeV 
(see Figs. 3, 10, 11), for tan/3 = 4 (left panel). Similar observations can be made for the right 
panel, in which case the maximal value of p% ~ 0.009% occurs at p = 500GeV (see Fig. 11), 
and like the former case, it decreases relatively with the increasing mass again. For instance, 
p 3 ~ 0.005%, at p, = 800GeV (see Figs. 3, 10, 11). Finally, it reaches far below p 3 ~ 0.001% for 
m h3 = 127GeV. 

As both panels of the figure shows that lighter the Higgs boson (/13) larger its CP-odd 
composition. Therefore, as the p^-m^ dependence (which corresponds to Q = 1000 GeV) 
suggests, any possible increase in the lower experimental bound of the lightest Higgs mass will 
imply reduced CP-odd composition, being in a reasonable agreement with the results of (3{| . 

Before concluding, we would like to note that we have concentrated on a particular scale 
Q = 1000 GeV, for which case the results are very similar to what we have found in the 
scale-independent case In fact, in both low and high tan/3 regimes, Q = 1000 GeV, is an 
approximately intersectional value of the second portion of the parameter space ( 925 < Q < 1200 
for tan/3 = 4, and 775 < Q < 1200 for tan/3 = 30), beyond which the p% — p dependence 
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changes. Moreover, if we consider the remaining interval in the second portion of the parameter 
space 1000 < Q < 1200, similar results can be obtained. Naturally, in the 1000 < Q < 1200 
interval, with the decrease in p, p% increases gradually, which is a characteristic behaviour of 
the second portion of the parameter space (see Fig. 7). However p 3 never exceeds ~ 0.02% even 
at Q = 1200 GeV. Therefore, one can conclude that there is a strong influence of Q on the 
Ps-p interdependence in the first portion of the parameter space (175 < Q < 925 for tan (3 = 4, 
and 175 < Q < 775 for tan/3 = 30). On the other hand, the scale dependence of the radiative 
corrections is sufficiently suppressed in the second portion of the parameter space, in particular, 
1000 J; Q < 1200, interval, for both low and high tan/3 regimes. 

4 Conclusion 

In this work, we have mainly concentrated on the percentage CP odd-compositions of the lightest 
Higgs boson, with the inclusion of the bottom-sbottom contributions, for all values of Q ranging 
from top mass to TeV scale. We would like to briefly summarize the results: 

1. A comparative glance at both low and high tan j3 regimes shows that as Q increases, the 
percentage CP-odd composition of the lightest Higgs (^3) decreases until Q = 925GeV 
for tan j3 = 4, and Q = 775GeV for tan (3 = 30. Beyond these points, p% gradually 
starts to increase in the admitted range of Q. However, it never exceeds 0.17% for all 
values of Q changing from 175 GeV to 1200 GeV. As compared to its percentage CP-even 
compositions this value is so small to cause observable effects. Therefore, lightest Higgs 
remains essentially CP even. 

2. The present analysis differs from that of [p9f |, in the sense that the latter does not depend on 
Q explicitely since the D-term as well as the bottom-sbottom contributions are not taken 
into the consideration. However, in the appropriate limit, namely for Q = lOOOGeV, a 



comparision between the results is possible, and a reasonable agreement with 3S] can be 
found. Moreover, similar observations can be made for the second portion of the parameter 
space, particularly in the 1000 < Q < 1200 interval, which is an allowed intersectional 
region in both low and high tan (3 regimes. Therefore, one can conclude that there is a 
strong influence of Q on the p% — p interdependence in the first portion of the parameter 
space. On the other hand, the scale dependence of the radiative corrections is sufficiently 
suppressed for the second portion of the parameter space, particularly in the 1000 ^ Q Z; 
1200 interval. However, consideration of various renormalization scales, all being around 
the weak scale, lead us to a wealth of CP violation opportunities. 
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3. The underlying model provides a quite restricted parameter space due to the naturalness 
requirements, as well as a simultaneous solution both to the the strong CP and //-problems. 

M. B would like to thank the Scientific and Technical Research Council of Turkey (TUBITAK) 
for partial support under the project, No:TBAG2002(100T108). 
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